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Abstract 

The discrete nonlinear Schrodinger equations of n sites are studied 
with periodic boundary conditions. These equations have n branches 
of standing waves that bifurcate from zero. Traveling waves appear as 
a symmetry-breaking from the standing waves for different amplitudes. 

The bifurcation is proved using the global Rabinowitz alternative in 
subspaces of symmetric functions. Applications to the Schrodinger 
and Saturable lattices are presented. 
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1 Introduction 

The discrete nonlinear Schrodinger equation (DNLS) appears in the study of 
optical waveguide arrays and Bose-Einstein condensates trapped in optical 
lattices [lO]. In this paper, we consider a general lattice of n sites described 
by the equations 


iqj = V'{\qjf)qj + {qj+i - qj) + {qj-i - qj) , ( 1 ) 
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where the sites, qj{t) G C for j = satisfy the periodic boundary 

conditions qj = qj+n- The equations include the Schrodinger lattice, V{x) = 
cx'^12, and the Saturable lattice, V{x) = cln(l +x), where c represents the 
strength of the linear coupling after rescaling. 

Equations ([T]) have the explicit solutions 

Pvr 

qj{t) = with C = —, (2) 


for m = 1,..., n, where the frequency a; is a function of the amplitude a G M"*" 
given in (j6]). These solutions are relative equilibria and are known as standing 
waves in the sense that their norms are stationary in time. 

We prove existence of periodic solutions using a non-abelian group that 
acts by permutating the oscillators and shifting and reflecting phase and time 
(see Definition [T]). In fact, solutions ([2]) appear from symmetry-breaking of 
the trivial solution; from them, we prove a secondary symmetry-breaking of 
periodic solutions (Theorem [5]). 

Main result. Let m G [0,n/2] flN with m ^ njA., for each k G [l,?7,/2] flN 
such that 

0fc(a) G (-CX), l)\{7i : j e [1, n/2] n M}, (3) 

the relative equilibrium (|2]) has two global bifurcations of solutions of the 
form 

q.{t) = [a + x {ut ± jkQ ), (4) 

where 


0fc(a) = 


2 cos mC sin^ 


and 


2kC 


7 fc = 1 — cot tan m(. 


Each branch is a global continuum in the space of 27r-periodic functions x 
and frequencies u emanating from (0, given in flT7|) . 

These solutions are traveling waves in the sense that their norms satisfy 


\qj\ (f) = a + r{ut±jkC), 

where r{t) is real 27r-periodic; these solutions are known as traveling or mov¬ 
ing breathers when they are localized. The existence of localized traveling 
waves in infinite lattices is proved in [IH] (see also Chapter 16 and 21 in [TO] 
and [3]). 
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In Theorem [7] we prove that solutions ([2]) are stable if the conditions ([2]) 
hold for k = 1,n — 1. In HD], solutions given by ([2]) are called plane waves, 
the nonlinear dispersion relation (6.7) is equivalent to (| 6 ]) and the modulation 
stability ( 6 . 8 ) has stable directions precisely when (|3]) holds. 

In [7], the authors hnd a bifurcation of relative equilibria for m = 1 and 
amplitudes 0fc(a) = 7 ^. This bifurcation exists due to an eigenvalue of the 
linearization that crosses zero. This phenomena occurs for any m, and we 
should expect a bifurcation of standing waves appearing from the amplitudes 

0fc(a) = 7A:. 

For (j)k{ci) > 1, there are two eigenvalues colliding on the imaginary 
axis and detaching into the complex plane. This phenomenon may trig¬ 
ger a Hamiltonian-Hopf bifurcation where isolated traveling waves persist 
for 0fc(a) > 1. Indeed, a bifurcation of this kind is described in [9] for the 
trimer n = 3. 

The authors prove in [H [ 6 l |5] bifurcation of periodic solutions for bodies, 
vortices and Schrodinger sites, for m = 1. In the body and vortex problems, 
the coupling is homogeneous and invariant under all permutations, i.e. the 
stability and bifurcation properties of ([ 2 ]) are independent of a and m; but this 
is not the case for Schrodinger sites. We complete the analysis of bifurcation 
and stability for all m’s. 

References glElls] use a global Lyapunov-Schmidt reduction and a topo¬ 
logical degree for G-equivariant maps that are orthogonal to the generators 
(see III Eli). In this paper we present a direct and self-contained approach, 
requiring non-abelian group actions and the global Rabinowitz alternative 

ra- 

In Section 2, we dehne the equivariant properties of the bifurcation oper¬ 
ator, we present a reduction to a hnite number of Fourier components, and 
we hnd the spectra. In Section 3.1 we analyze the spectra for each irreducible 
representation and we prove the bifurcation result. In Section 3.2 we study 
the stability. In Section 4, we apply the results to the focusing and defocusing 
Schrodinger and Saturable lattices. 

2 Setting the problem 

Equations ([T]) in rotating coordinates, qj(t) = e^‘^^Uj(t), are 

iiij - uuj = V\\ujf)uj + {uj+i - Uj) + {uj-i - Uj ). (5) 
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The values aj = satisfy aj = aj+n and 

(oj+i — 2aj + = (—4sin^ mC/2)aj. 

Then Uj(t) = aj is an equilibrium and (|2]) is a solution of ([I]), when 

(n = 4sin^mC/2 — Id'(a^). (6) 

In real coordinates, Uj G equations (jS]) are 

Jiij = ojUj + V'(\uj\^)uj + (uj+i — 2uj + Mj-i), 
where J is the symplectic matrix 

J=('; -;) andfl=(J _°i), (7) 

Let u = (wi, ...,Un) and J' = diag{J,J). The system of equations in 
vectorial form is 

Jii = WHiu), 

where H is the Hamiltonian 
1 ” 

H = - ^{V{\uj\^) + uj\uj\^ - \uj+i-Uj\^] . (8) 

i=i 

In this setting, the relative equilibrium is 

= («!, ■■■iCin) with aj = (9) 

Using the change of variables u{t) = + xi^ut), 27r/z/-periodic solutions 

of the Hamiltonian system correspond to zeros of the operator 

/(x; u) = Jx- + x) : x M+ ^ 

Below we will prove global bifurcation of periodic solutions from the set of 
trivial solutions (0, u) for u G 

Definition 1 Letm be as in The linear action of the non-abelian group 

G={ZnX S^) U K{Zn X S^) 
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in is given by the homomorphism of groups p : G ^ GL{Ll^) 

generated by 


p(C, ^)xj{t) = e + if), 

p{K)Xj{t) = RXn-j{-t), 

where ( = 2n/n generates in = [0,27r) and k is the reflection that 
generates two copies ofl^n x 5'^ in G. 

Since the Hamiltonian H is gauge invariant and autonomous, the map / 
is G-equivariant. Moreover, all the elements of G leave the equilibrium 
hxed, and then the isotropy group of is Ga„, = G. 

2.1 Irreducible representations 

In this section we find the irreducible representations of the action of G. We 
dehne the isomorphisms —)■ 14 as 

TkZ = ^ 

where 14 is the image of for k G {1, 

A function x{t) G L 2 ^(R^") is expressed in the Fourier components as 
X = Fourier component Xi G as Xi = Ylk=i'^kXk,i 

with Xk,i G C^. Therefore, a function can be expressed in the ortho normal 
coordinates TkXk^ie^’’^ as 

x(t) = ^ TkXk,ie"^\ = {1, (11) 

(k,l)£Z„xZ 

Let us denote the j-th component of TkXk,i G by the two-dimensional 
vector G C^. With this notation, the j-th component of 

p{C)TkXk,i is 

then 

p{C)TkXk,i = e"^^TkXk,i. 

Similarly, 

p{K)TkXk,i = TkRxk,i. 
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Therefore, the subspaces of similar irreducible representations of the space 
by the action of G are 

= {TkXk,iG^^ '■ Xk,i e C^}. 

In the components Xk,i, the action of G is 

p{C)xk,i = G^‘^Xk,i, p{p)xk,i = G^‘^Xk,i, p{K)xk,i = Rxk,i. ( 12 ) 

2.2 Linearization 

Since Vk^i are the subspaces of similar irreducible representations of G and 

/'(O) = Jdt - 

is G-equivariant, by Schur’s lemma, the linearization f'{0) is block diagonal in 
the components 14^;. The diagonal decomposition can be obtained explicitly 
using the following proposition, 

f’{0)x= Tk{UJ - Bk)xk,iG^\ 

(fc,z)ez„xz 

Proposition 2 Let ak and (3k he 

ttfc = dcosm^sin^ A;^/2, = 2 sinmC sin(13) 

the Hessian satisfy that 

D^H{a^)TkZ = TkBkZ, 
where Bk is the 2x2 matrix 

Bk = diag(2aV"(a^) - Ofc, -ak) + (14) 

Proof. We express D‘^H{sLm) in 2 x 2 blocks Aij as 

= (A,,)i:,=i. 

Since the coupling in the DNLS equations happens only between adjacent 
sites, then Aij = / for |i — jj = 1 and Aij = 0 for \i — j\ > 1, modulus n. 
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Using ttj = ae-^™'^'^ei, we have 

^DV(a„) = U'(a 2 )/ + 2 aV"(a 2 )e(^’”«)^eie[e-(^”^^)^, 

where eiej = diag{l,0). Since U'(a^) + u — 2 = —2cosmC, we conclude 

Aj^j = -2(cosmC)/ + 2aV"(a2)e(^'”^^)^eiefe-(^™^)^. 

Given that the j-th component of TkZ E is E C^, 

where 2 ; E C^, the j-th component of D^H{a)TkZ is 

^ ^{ikI+mJ)C _j_ g-(ifc/+mJ)C^'\ ^j{ikI+mJ)C, ^ ^ ^j(ikI+mJ)C 

where 

Bk = -2(cosmC)/ + 2a^V"{a^)eieJ + + e-h^-f+”^UC, 

From the equalities 

^iiki+mj)c ^ ^-iiki+mj)( ^ (2cosA;CcosmC)/ + (2sinA;CsinmC)f J, 

and 

—2cosmC(l — cosk() = —Acosm({smk(/2Y = —at, 
we conclude that 

Bk = -aki + /3k{iJ) + 2a^V"{a‘^)diag{l, 0 ). 


3 Main Results 


The symmetries permit us to assume, without loss of generality, that m E 
[0,n/2] n N. Throughout this section we also assume m 7 ^ n/4. 

In these cases the sign of is well dehned for k = 1,..., n — 1, 


sgn{ak) = 


1 if m e [0, n/4) 

— 1 if m e (n/4, n/2] 


1 
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and we can define 

Ma) = —V"(a^), 7 t = l-f*y. (15) 

Oik \OikJ 

Therefore, the matrix D‘^H{aim) is block diagonal with blocks 

Bn = diag{2a^V"{a^), 0), = akdiag{4>k - 1, -1) + /3k{iJ), 

ioT k = 1, ..,n — 1 . 

3.1 Bifurcation theorem 

We consider bifurcation in the hxed point space of the isotropy group Dn 
generated by (C, —kQ and k. Solutions in the hxed point space of Dn have 
symmetries 

Xjit) = P(C, -K)xj{t) = - kC). ( 16 ) 

That is Xj(t) = e^^^^Xnit + jkQ and, by the action of k, 

Xn{t) = p{K)Xn = RXn{-t). 

The component Xk,i G is hxed by k G -D„ if Xk,i G M x iM. 

Lemma 3 The matrix iJB^ in the subspaee M x iR has eigenvalues 

u^ = I3k± ^Jal{l-(j)k). (17) 

Moreover, 

(a) If k E [l,n — 1] nN and 0fc(a) G {—oo,'yk), then is positive. 

(b) If k E [l,n/2) fl N and (fkici) £ {lk, 1)? then irf and are positive. 

Proof. Let L = diag(l,i). The eigenvalues of iJB^ in the subspace R x iR 
are the eigenvalues of the real matrix 

L-\i.IBk)L={ • 

^ ’ \oik {(fk - 1 ) (Ik J 

The eigenvalues u of this matrix are the zeros of 

dkii^) = (^^ - I3kf - Oil{l - (fk)- 
8 


( 18 ) 



The function dk{0) has two real solutions if and only if 0^ G (—oo,l). 
Since dk{i>) is a polynomial of order at inhnity, then is positive and 
is negative if dk{0) is negative. This is the case if 0^ < 1 — {(dk/akY = Ik- 
For 0 G (7 a:, 1), the function dA:(0) is positive, and the values have the 
same sign of (3k- We conclude this result from the fact that 0*; > 0 for 
k G [1, n/2) n N. ■ 

Definition 4 We say that the amplitude a is non-degenerate ifV"{o?) ^ 0, 
0fc(a) Y Ik for k = — 1. We say that the frequency vf for k = 

1,n — 1 is non-resonant if uf ^ luf for integers I > 1 and j ^ k. 

The matrix Bn restricted to the subspace Xnp = p{n)xn,o = Rxn,o has the 
simple eigenvalue 2a^l/"(a^). For k = the blocks Bk have determi¬ 

nants 0| — al{l — 0fc). Therefore, the non-degeneracy property of a assures 
that D'^Hfsinf) has no zero-eigenvalues in Fix(£)„). 

Theorem 5 Assume that a is non-degenerate. If is non-resonant, then 

(a) For each fc G [1, n — 1] fl M such that 0fc(a) G (—oo, jk), there is a global 

bifurcation from (0, uf ) in the space 

{x G : Xjft) = e^'^^-^Xn{t +jkC,),Xn{t) = Rxn{-t)} X R+ (19) 

(b) Fork G [l,n/2]nN such thatcfkio) G (y/c, 1), there are global bifurcations 

from (0, z/^) and (0, z/0) in the space given by 

Proof. Let K : —)■ be the operator dehned in the Fourier basis 

XE j}f 

'T* , Q C 

Kx = Xo + {lij)~^xie^^^. 

ZGZ\{0} 

Since K : —>■ is compact and 

f{x) = Jx — i'~^D‘^H{aLm)x + 0{ |a;|^), (20) 

then 

Kf{x) =x- T{iy)x + g{x) : (21) 

where T(z/) is the compact linear operator 

Tx={I + W^D^H{si„,))xo + {yl)-\jD^H{sim)xiF^\ 

zgz\{o} 
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and g{x) = 0{\x\^^i ) is & nonlinear compact operator. 

Since / and K are G-eqnivariant, the operator Kf is G-eqnivariant. Then 
Kf is well dehned in the space Fix(Zi)„) given in flT^ . The global bifnrcation 
follows from Theorem 3.4.1 in [T3] if T[v) has a simple eigenvalne crossing 1 
in Fix(T)„). 

For Z = 0, by the non-degeneracy property of a, the matrix I+v~^D^H (a^) 
has no eigenvalnes eqnal to 1 in Fix(Z)„). For I > 2, dne to the non-resonance 
property of the matrix H(abm) has no eigenvalnes eqnal to 

1. For / = 1, we have that has an eigenvalne crossing 1 when 

z/ crosses z/^, corresponding to the block iJBk. Moreover, this eigenvalnes is 
simple when iJB^ is restricted to M x zM. We conclnde that T{v) in Fix(Zl„) 
has a simple eigenvalne crossing 1 when u crosses vf. ■ 

Remark 6 In the case of 1 \ I resonances, for I > 2, the previous 

theorem gives the existence of the bifurcation with the biggest frequency uf. 
In the case of 1 : 1 resonances, cannot prove existence of 

bifurcation because there is a double eigenvalue ofT{y) crossing 1. In this 
case, a Hamiltonian-Hopf bifurcation may appear at (fkia) = 1- This is 
described in Theorem 11.5.1 of fT^ . where two isolated solutions persist for 

0fe(a) > 1. 

3.2 Stability Analysis 

Let 

_ f sgn{V"{a‘^)) if m e [l,n/4) 

\ —sgn{V"{a‘^)) ii m e {n/4,n/2] 

Theorem 7 If am < 0, or am > 0 and 0i(a) < 1, then the relative equilib¬ 
rium is linearly stable. 

Proof. The Hamiltonian eqnation (|8]) is linearly stable at the equilibrinm 
a^ if JD‘^H{aim) has only pairs of pnrely conjngated imaginary eigenvalnes 
±zz/. Since J'D‘^H{am) has a pair of zero eigenvalnes that corresponds to the 
gange symmetry, the system is linearly stable if J’D'^H^a.m) has rz, — 1 pairs 
of pnrely imaginary eigenvalnes. 

The sign of (fk does not depend on fc G {1,..., n—1} and is eqnal to am- For 
am negative, we have 0fc(a) < 1. For am positive, nsing the fact that cos kC. 
is increasing for k E [1, rz,/2] fl N, then (fk is decreasing for k = 1,..., rz,/2 and 
0fc(a) < 0i(a) < 1. In both cases, it is not difficnlt to see that J'D‘^H{sLm) 
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has eigenvalues iu^ for fc = 1,n — 1. Since we conclude that 

J'D‘^H{airn) has n — 1 pairs of purely imaginary eigenvalues. ■ 


4 Applications 

In complex coordinates, 

q.{t) = (a + + jkQ). 

These solutions are discrete traveling waves in the sense that the norms 
satisfy 

\qj\ (t) =a + \xn\ (ut + jkC). 

For example, if k divides n, the traveling wave has k identical waves with 
wavelength equal to n/k sites. 

In the following sections, we present applications to the Schrodinger and 
Saturable lattice. 

4.1 Schrodinger lattice 

The cubic Schrodinger potential is V{x) = ca:^/2, where c > 0 corresponds 
to the focusing case, and c < 0 to the defocusing case. Standing waves given 
by ([2]) exist for 

a; = 4 sin^ nr^/2 — ca?. 

In the focusing case, the potential satishes 

V''{a?) = c > 0 and a = 1. 

If m G [0, n/4), then ctfc > 0 and 0fc(a) < 1 for < ak/2c. If m G (n/4, n/2], 
then ttfc < 0 and 0fc(a) < 1 for all a. 

Therefore, the following result holds. 

Proposition 8 In the focusing Schrodinger lattice (c > 0), the equilibrium 
S-m is linearly stable when m G (n/4, n/2], or m & [0,n/4) and a < y/oiplc. 
Moreover, the equilibrium has two global bifurcations of traveling waves 
for each k G [1, n/2] flM z/m G (n/4, n/2], ormE [0, n/4) and a < yjak/2c. 
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In the defocusing case, the potential satishes 

= c < 0 and a = —1. 

For m e [0,n/4), 0fc(a) < 1 . For m G (n/4,n/2], if then 

0fc(a) < 1 . 

Proposition 9 In the defocusing Schrodinger lattice (c <Q), the eguilibrium 
Am is linearly stable when m G [0,n/4), or m & (n/4, n/2] and a < a/q;i/2c. 
Moreover, the eguilibrium Am has two global bifurcations of traveling waves 
for each k G [1, n/2] flN i/m G [0, n/4), ormE (n/4, n/2] and a < 

4.2 Saturable lattice 

The Saturable potential is given by V{x) = cln(l +a;) with c > 0. It is clear 
that 

V'\a^) = —c(l + and a = —1. 

If m G [0,n/4), then 0fc(a) < 1 for all a, while if m G (n/4,n/2], then 
0fc(a) < 1 for (a + a“^) < —ak/‘2.c. 

Proposition 10 The eguilibrium Am in the Saturable lattice is linearly sta¬ 
ble when m G [0,n/4), or m E (n/4,n/2] and (a + a“^) ^ < \J—ail2c. 
Moreover, the eguilibrium Am has two global bifurcations of traveling waves 
for each k E [1, n/2] fl N if m E [0, n/4), or m E (n/4, n/2] and 

(a + a“^) ^ < a/—Q; fc/2c. 

Acknowledgements. C. Garcia is grateful to P. Panayotaros, M. Tejada- 
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improved the presentation of this manuscript. 
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